Abstract Sounding data are first designed to be applicable to wide-angle refraction/reflection, and the 1D velocity structure linear inversion of iterative calculation is conducted by using travel time residual as the threshold. Then LIAE1D is programmed accordingly with the following characteristics: 1) the establishment of an adaptive preliminary model| 2) the solution of the linear equation by means of a broader inverse moment matrix. The method and program are later proved to be feasible and efficient in numerical simulation tests as well as in the application of the Huian-Datian-Ninghua profile of Fujian artificial earthquake sounding wide-angle refraction/reflection. Thus, these results could provide a suitable preliminary model for future 2D inversion.
